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Abstract 
Several years back Angliker et al. [Chem. Phys. Lett. 87, 208 (1982)] predicted nonacene 
to be the first linear acene with the triplet state l 2, B 2u as the ground state, instead of the 
singlet l x A g state. However, contrary to that prediction, in a recent experimental work 
Tonshoff and Bettinger [Angew. Chem. Int. Ed. 49, 4125 (2010)] demonstrated that 
nonacene has a singlet ground state. Motivated by this experimental finding, we decided 
to perform a systematic theoretical investigation of the triplet states of polyacenes, with 
an emphasis on singlet-triplet splitting, starting from naphthalene, all the way up to de- 
cacene. Methodology adopted in our work is based upon Pariser-Parr-Pople model (PPP) 
Hamiltonian, along with large-scale multi-reference singles-doubles configuration interaction 
(MRSDCI) approach. Our results predict that the singlet-triplet gap decreases with acene 
size and octacene onwards, it reaches saturation. Thus no singlet-triplet crossover takes 
place in longer acenes and the ground state is of singlet multiplicity. We also analyze the 
nature of many-particle wavefunction of the correlated singlet ground state and find that the 
longer acenes exhibit tendency towards a open-shell singlet diradical ground state. More- 
over, when we compare the experimental absorption spectra of octacene and nonacene with 
their calculated singlet and triplet absorption spectra, we observe very good agreement for 
_L the singlet case. Hence, the optical absorption results also confirm the singlet nature of the 

ground state for longer acenes. Calculated triplet absorption spectra of acenes predict two 
well separated intense long-axis polarized absorptions, as against one such peak observed 
for the singlet case. This is in important prediction regarding the triplet optics of acenes, 
, which can be tested in future experiments on oriented samples. 
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I. INTRODUCTION 

Polyacenes, which can be seen as linearly fused benzene rings, are known for their well 
defined structures, and crystalline forms. 1 2 Because of their small band gaps and high charge- 
carrier mobilities, they find potential applications in novel opto-electronic devices such as 
light-emitting diodes, and field effect transistors etc., which make them experimentally and 
theoretically a very important class of materials.^ In spite of a long tradition of research,^ 
the field of acenes has experienced a resurgence of interest in recent years because they 
are also perceived as the building blocks for organic electronic materials such as graphene 
nanoribbonsP^ 

Although pentacene, has excellent optical and transport properties, however, it is con- 
ceivable that the longer acenes could have even more attractive properties, with possible 
applications in the field of nanotechnologyP As the size of the longer acenes approaches the 
nanometer scale, their reactivity also increases, and, therefore, it has been difficult to syn- 
thesize them from heptacene onwards.^ Recently, many efforts have been made to synthesize 
the longer acenes, e.g., heptacene and its functionalized derivatives have been synthesized by 
several workersP^ By using cryogenic matrix-isolation technique, and a protection group 
strategy, octacene and nonacene have been synthesized by Tonshoff and Bettinger^, while 
Kaur et al™ prepared functionalized nonacene. 

While all the known oligoacenes ranging from naphthalene to hexacene have a singlet 
ground statep^ some years back, based on an extrapolation of the experimental singlet- 
triplet gap of the acenes up to pentacene, Angliker et alP^* predicted nonacene to be the 
first linear acene with a triplet state (l 3 i?2«) as the ground state, instead of the singlet one 
(l 1 ^). This was an interesting prediction because, if true, it could open the possibilities of 
magnetic applications of longer acenes. However, contrary to that prediction, based upon 
their recent optical absorption experiments, Tonshoff and Bettinger^ demonstrated that 
nonacene has a singlet ground state. Motivated by this experimental finding, we decided to 
perform a systematic theoretical investigation of the electronic structure of the ground and 
low-lying excited states of longer acenes, with an emphasis on the singlet-triplet splitting, 
and their optical properties. In order to realize the possible potential of longer acenes in nan- 
otechnology, a deep theoretical understanding of their electronic structure is very important. 
For our calculations, we adopt a methodology based upon Pariser-Parr-Pople model (PPP) 
Hamiltonian, along with large-scale multi-reference singles-doubles configuration interaction 
(MRSDCI) approach. We present results on the singlet and triplet low-lying excited states 
of octacene, nonacene, and decacene, along with their optical absorption spectra. The cal- 
culated spectra are compared with the experimental ones for octacene and nonacene, and 
excellent agreement is found between our singlet optical absorption spectra, and the exper- 
imental ones. This, along with our results on singlet-triplet splitting, leads us to conclude 
that the acenes in question possess singlet ground states. We also make important predic- 
tions about triplet absorption in long acenes, which can be tested in future experiments. 

The remainder of this paper is organized as follows. In Section [IT] we briefly review the 



theoretical methodology adopted in this work. In Section III we present and discuss our 
results of: (a) the singlet-triplet splitting of oligoacenes of increasing lengths, (b) singlet 
linear absorption spectra, (c) triplet optical absorption spectra, and (d) the comparison of 
the experimental absorption spectra of octacene and nonacene, with their calculated spectra. 



Finally, in Section IV we summarize our conclusions, and present possible directions for the 



future research work. 



II. THEORY 

The schematic structures of higher polyacenes studied in this work are shown in [TJ The 
molecule is assumed to lie in the xy-plane, with the conjugation direction taken to be along 
the x-axis. The carbon-carbon bond length has been fixed at 1.4 A , and all bond angles 
have been taken to be 120°. The reason of choosing this symmetric geometry, against various 
other possibilities has already been discussed in our earlier paperP^^ It can be noted that 
these structures can also be seen as two polyene chains of suitable lengths, coupled together 
along the y-direction. 
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Figure 1. Schematic drawings of long acenes (octacene to decacene, from top to bottom) considered 
in this work. The x axis is assumed along the conjugation direction, while the y axis is perpendicular 
to it, in the plane of the figure. 



The correlated calculations are performed using the PPP model Hamiltonian, which can 

be written as 

H = H Cl + H c , 2 + H Cl c 2 + H ee , 

where Hc 1 and Hq 2 are the one-electron Hamiltonians for the carbon atoms located on the 
upper and the lower polyene like chains, respectively. Hc 1 c 2 * s the one-electron hopping 
between the two chains, and H ee depicts the electron-electron repulsion. The individual 
terms can now be written as, 

{k,k') 
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In the equation above, k, k' are carbon atoms on the upper polyene chain, /i, v are carbon 
atoms located on the lower polyene chain, while % and j represent all the atoms of the 
oligomer. Symbol (...) implies nearest neighbors, and B it j = ^2 ff (c\ a c^ a + h.c), where 

c ia (cj )CT ) denotes the creation (annihilation) operator for a tt orbital of spin a, located on 
the ith carbon atom. Matrix elements t , and t± depict one-electron hops. As far as the 
values of the hopping matrix elements are concerned, we took to = 2.4 eV for all nearest 
neighbor hopping with t± = to, consistent with the undimerized ground state for polyacene 
argued by Raghu et air® 

The Coulomb interactions are parametrized according to the Ohno relationship p*l 

V l , j = U/K lJ (l + 0.QH7Rl J ) 1/2 , 

where, Kij depicts the dielectric constant of the system which can simulate the effects of 
screening, U is the on-site repulsion term, and Rij is the distance in A between the ith and 
the jth carbon atoms. In the present work, we have performed calculations using "standard 
parameters'^ with U = 11.13 eV and K iy j = 1.0, as well as "screened parameters'^ with 
U = 8.0 eV and n it j = 2.0 (i ^ j) and K iti = 1.0. The screened parameters employed here 
were devised by Chandross and Mazumdar 20 with the aim of accounting for the inter chain 
screening effects in phenylene based polymers. Our motive behind using these parameters 
is precisely the same. 

The starting point of the correlated calculations for the molecules is the Restricted 
Hartree-Fock (RHF) calculations, employing the PPP Hamiltonian, using a code devel- 
oped in our group.^ All the resultant HF molecular orbitals are treated as active orbitals. 
The many-body effects beyond RHF are computed using the Multi-Reference Singles Dou- 
bles Configuration Interaction (MRSDCI) methocP^"^ in the following manner. After RHF 
calculations of the ground state l l A g , the MRSDCI calculation of the ground state, l^Ag 
and excited states of symmetries l B 2u and l B 3u are performed by taking the lowest energy 
configuration of the corresponding symmetry (lB 2u is H — > L and lB 3u is H — )■ L + 1 and 
H — 1 — > L , where H and L corresponds to Highest Occupied Molecular Orbital (HOMO) 
and Lowest Unoccupied Molecular Orbital (LUMO)). A similar approach is adopted for 
the triplet states as well, except that one has to ensure that the multiplicity of the wave 
functions is triplet. From the CI calculations, we obtain the eigenfunctions and eigenval- 
ues corresponding to the correlated ground and excited states of the examined molecules. 
Using these eigenfunctions, the dipole matrix elements between the ground state and vari- 
ous excited states are computed. For the triplet states, the matrix elements are calculated 
with respect to the lowest triplet state l 3 B 2u - These dipole matrix elements, along with 
the energies of the excited states are, in turn, utilized to calculate linear (or triplet) optical 
absorption spectrum. Important excited states contributing to various peaks of the spec- 
trum are identified, and a new set of MRSDCI calculations are performed with an increased 
number of reference configurations contributing both to the ground state, and the excited 
states, leading to a new absorption spectrum. This procedure is iterated until the computed 
spectrum converges satisfactorily. 

III. RESULTS AND DISCUSSIONS 

In this section we present the results of our MRSDCI calculations on polyacenes ranging 
from naphthalene to decacene examining their singlet-triplet gaps, with the aim of deter- 
mining the spin-multiplicity of their ground state. Furthermore, we also present calculations 



on the optical absorption of the long acenes, namely, octacene to decacene from their lowest 
triplet states (l 3 B 2u ), and the lowest-singlet state (l 1 ^), and compare the results with the 
experimental ones, where available. 



A. Singlet-Triplet Splitting 

The first set of calculations which we performed to explore singlet-triplet (l 1 Ag-l 3 B 2u ) 
splitting in oligoacenes were at the RHF level for the singlet states, and the UHF (Un- 
restricted Hatree Fock) level for the triplet states. For an oligomer containing 2m elec- 
trons, in the UHF calculations m + 1 electrons were assigned the majority spin, and the 
remaining m — 1 electrons the minority spin to represent the triplet ground state. The 
results of these calculations for oligoacenes up to acene-15 (acene-n denotes an oligomer 
with n benzene rings) are presented in III Examining the values of the singlet-triplet split- 
ting AE = E(l 3 B 2u ) — Eil^Ag) for various oligomers, it is obvious that: (a) in standard- 
parameter based calculations, acene-10 onwards triplet state become lower in energy as 
compared to the singlet state, while, (b) with screened parameters, this crossover takes 
place at smaller lengths, starting with acene-7. Therefore, these results appear to confirm 
the essential prediction of Angliker et a/P^ that the singlet-triplet crossover will take place 
in oligoacenes with the increasing lengths, although, they differ in details in that Angliker 
et al. 15 predicted the crossover from acene-9 onwards. However, the shortcomings of these 
calculations are that they have been performed at the mean-field level, and, therefore, do 
not include the electron-correlation effects. Although, the calculations of Angliker et alP^ 
were performed using the CI approach, but they employed the singles-CI (SCI) method, 
which accounts for the electron-correlation effects in a very limited manner. Furthermore, 
another calculation, performed by Houk et al. p^ also predicted the singlet-triplet crossover 
in oligoacenes. However, that calculation was performed using an ab initio DFT based 
approach,^ and, therefore, did not account for the electron-correlation effects. Based upon 
past calculations, 27 it is a well-known fact that in order to be able to predict the correct 
excited state orderings in conjugated polymers, it is very important to account for the 
electron-correlation effects in an accurate manner. Therefore, we decided to go beyond the 
HF approach, and performed large-scale CI calculations to explore the singlet-triplet order- 
ing in polyacenes. Before we present and discuss our results, we would like to give a flavor to 
the reader as to the size of the CI calculations performed. [IT] presents the number of reference 
states used in the MRSDCI calculations, and the total size of the resultant CI matrix. The 
largest calculation consisted of more than one million configurations, from which it is obvious 
that our calculations were large scale, and, thus, they account for the electron-correlation 
effects in an adequate manner. The calculated values of the singlet-triplet gap, and their 
comparison with the experimental and other theoretical calculations, are presented in |IH[ for 
the oligoacenes ranging from naphthalene to decacene. The variation of the singlet-triplet 
gap, as a function of the oligomer length, is presented in [2j It is obvious both from the 
table, and the figure, that the excitation energy of l 3 B 2u state decreases as the oligomer 
length increases. Nevertheless, even for decacene the singlet-triplet gap is nonvanishing, and 
appears to saturate as a function of the increasing chain length. As is obvious from the [2] 
that similar trends were also observed in the works of other authors, except the DFT-based 
UB3LYP calculation performed by Bendikov et alP^ which reports much smaller singlet- 
triplet gaps. We believe that these smaller gaps obtained in the work of Bendikov et aZ.pS 
could be attributed to the well-known tendency of DFT to underestimate the energy gaps. 



Therefore, based upon the fact that for decacene our screened parameter based calculations 
predict a singlet-triplet gap ps 0.5 eV, we conclude that even for longer acenes, the singlet 
state l 1 ^ will be the ground state, and thus no singlet-triplet crossover of the kind predicted 
by Angliker et al.p^ occurs as per our calculations. As far as comparison with other au- 
thors is concerned, our screened parameter results are in good quantitative agreement with 
DMRG/STO-3GP results, and good qualitative agreement with the DMRG/DZ resultsP 
Recently, Hajgato et a/P^ performed first principles coupled cluster CCSD(T) calculations on 
the singlet-triplet gaps of acenes ranging from octacene to undecacene (acene-11), and their 
reported value of 0.58 eV for octacene is in perfect agreement with our screened parameter 
results. However, for nonacene and decacene their reported values 0.46 eV and 0.35 eV^, 
respectively, are smaller than our screened parameter values of 0.56 eV and 0.54 eV. Regard- 
ing the singlet-triplet gap in the polyacene limit (n — > oo), Raghu et al. pausing the DMRG 
method and the PPP model, predicted it to be 0.53 eV, while Gao et al., 32 using a spin 
Hamiltonian, estimated it to be 0.446 eV, both of which are close to our screened-parameter 
value of 0.5 eV computed for decacene. 



Bendikov et al. 28 noted that the restricted singlet density functional ground state of higher 
acenes would become unstable due to its open-shell nature, with two unpaired electrons (a 
singlet diradical 33 ) for acenes longer than hexacene. Based upon ab initio DMRG calcula- 
tions on acenes in the range n = 2-12, Chan and coworkers found that the ground state 
wave functions for longer acenes were of the type of polyradical singletsP^ In an another 
DFT work, Jiang and DaP^ predicted the ground state of octacene and higher acenes to 
be antiferromagnetic (in other words, open-shell singlet), but not necessarily a diradical. 
Unlike the DFT calculations, in our approach, the many-body wave functions of the ground 
state and the excited states of the studied oligomers are available. Therefore, we decided 
to probe the nature of the ground state of longer acenes to ascertain whether, or not, they 
exhibit a polyradical character. The character of the many-body wave functions of the l^Ag 
ground state, obtained in our best CI calculations on oligomers ranging from naphthalene to 



decacene are presented in IV, for both the standard, and the screened parameters. From the 
results, it is obvious that the contribution of the configurations with two open-shell orbitals 
increases with respect to the HF reference state with the increasing oligomer length, while 
the contribution of the closed-shell HF state decreases. Therefore, based upon our results, 
we conclude that the longer acenes do exhibit a tendency towards a singlet diradical ground 
state. 
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n 


std 


scr 


S 


T 


AE 


S 


T 


AE 


2 


-22.73 


-20.02 


-2.70 


-19.19 


-17.09 


-2.10 


3 


-31.98 


-30.36 


-1.61 


-27.17 


-26.01 


-1.16 


4 


41.16 


-40.14 


-1.02 


-35.10 


-34.48 


-0.62 


5 


-50.32 


-49.70 


-0.62 


-43.02 


-42.72 


-0.29 


6 


-59.47 


-59.10 


-0.37 


-50.94 


-50.85 


-0.08 


7 


-68.62 


-68.41 


-0.21 


-58.85 


-58.89 


0.04 


8 


-77.76 


-77.67 


-0.10 


-66.76 


-66.89 


0.13 


9 


-86.91 


-86.88 


-0.02 


-74.67 


-74.86 


0.18 


10 


-96.05 


-96.07 


0.02 


-82.59 


-82.81 


0.22 


11 


-105.19 


-105.24 


0.05 


-90.50 


-90.74 


0.24 


12 


-114.34 


-114.41 


0.07 


-98.41 


-98.67 


0.25 


13 


-123.48 


-123.56 


0.08 


-106.32 


-106.59 


0.26 


14 


-132.62 


-132.71 


0.09 


-114.23 


-114.50 


0.27 


15 


-141.76 


-141.86 


0.10 


-122.15 


-122.42 


0.27 



Table I. RHF (singlet, S) and UHF (triplet, T), total energies along with their difference in energies 
(AE), in eV, of acenes n=l-15 computed using the standard (std), and the screened (scr) parameters. 



Table II. The number of reference configurations (N re t) and the total number of configurations 
{Ntotal) involved in the MRSDCI (or FCI or QCI, where indicated) calculations, for different sym- 
metry subspaces of various oligoacenes. 



n 


4, 


l B 2u 


3 B 2u 


N re f 


Ntotal 


Nref 


Ntotal 


N ref 


Ntotal 


2 


l a 


4936 a 


r 


4794a 


l a 


4816 a 


3 


l a 


623576 a 


i a 


618478 a 


l a 


620928 a 


4 


l b 


193538 a 


i b 


335325 a 


100 c 


323063 c 












86 d 


319005 d 


5 


l b 


1002597 6 


i b 


1707243° 


52 c ' d 


581702 c ' d 


6 


100 c 


1110147 c 


100 c 


1173212 c 


65 c 


1461526 c 




100 d 


1177189 d 


100 d 


1328252 d 


63 d 


1551590 d 


7 


35 c 


856788 c 


30° 


674925 c 


33 c 


1369624 c 




22 d 


QlttW d 


30 d 


850627 d 


29 d 


1300948 d 


8 


18 c 


768641 c 


14 c 


509119 c 


19 c 


1066355 c 




12 d 


540651 a! 


4 d 


145978 d 


14 d 


918645 d 


9 


13 c 


959737 c 


13 c 


769387 c 


18 c 


1626229 c 




12 d 


871397 a! 


3 d 


186651 d 


12 d 


1152071 a! 


10 


ll c 


1202681 c 


\T 


1192394 c 


15 c 


1735352 c 




U d 


1199887 d 


3 d 


270187 d 


10 d 


msise^ 



a FCI method with standard as well as screened parameters, 
b QCl method with standard as well as screened parameters, 
c using standard parameters, 



using screened parameters. 



Table III. For various oligomers, the comparison of the calculated singlet-triplet gaps (excitation 
energies of l 3 B 2u , with respect to the 1 A g ) with the experiments, and other theoretical results. 
The column std (scr) corresponds to the results obtained using the standard (screened) parameters 
in the PPP model Hamiltonian. 



n 


E(l 3 B 2u ) - E(l l A g ) (eV) 


std 


scr 


Expt. 


Other Theoretical 


2 


2.53 


2.11 


2.64^1 


2.71126] 2.73281 2.67291 


3 


1.73 


1.48 


lp8 -jBffl 


1.81126] 1.81128] L9 c|29] 


4 


1.25 


1.11 


1.27^ 


1.2(pl 1.16^1,1.511221 


5 


0.99 


0.93 


0.86 ± 0.03^ 


0.782S1, 0.76^1,1.16^1 


6 


0.87 


0.85 


0.54±0.05 15 


0.462S1, 0.4T«2SI,0.9ll22l 


7 


0.73 


0.69 


- 


024261 


8 


0.68 


0.58 


- 


0.0^61, 0.13^1,0.62^1,0.58^1 


9 


0.60 


0.56 


- 


0.07^1,0.50^1,0.46^1 


10 


0.57 


0.54 


- 


0.iP )0 .4P,0.3P 
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1 1 1 


i i i 
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■-■ Std. par. 
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Figure 2. (Color online) Singlet-triplet energy gap as a function of the acene lengths. 



Table IV. Magnitudes of the coefficients of the closed-shell Hartree Fock (HF), and the open shell 
(OS) doubly excited virtual configuration \H — > L; H— 1 — > L+l) to the singlet l l A g ground state 
CI wave functions of acene-n, obtained using the standard (std), and the screened (scr) parameters 
in the PPP model. In the preceding expression H stands for the highest occupied molecular orbital 
(HOMO), while L stands for lowest unoccupied molecular orbital (LUMO). 



n 


HF 


OS 


std 


scr 


std 


scr 


2 


0.94 


0.90 


0.14 


0.13 


3 


0.91 


0.86 


0.13 


0.11 


4 


0.88 


0.83 


0.11 


0.13 


5 


0.84 


0.79 


0.13 


0.14 


6 


0.81 


0.77 


0.15 


0.16 


7 


0.80 


0.76 


0.17 


0.19 


8 


0.79 


0.74 


0.18 


0.20 


9 


0.77 


0.73 


0.19 


0.21 


10 


0.76 


0.71 


0.20 


0.22 



B. Singlet Linear Optical Absorption calculations 



In an earlier work in our group, we had reported the calculations of linear optical ab- 
sorption spectra of oligoacenes ranging from naphthalene to heptaceneJ-^ In this work, we 
extend our calculations to longer acenes, and present the calculations of linear optical ab- 
sorption in octacene, nonacene and decacene. In |4j [5] and [6j we present the singlet linear 
absorption spectra of these oligoacenes from their ^- l A~ singlet state computed using the 
standard parameters and the screened parameters. As per dipole selection rules for the D^h 
symmetry, allowed one-photon transitions from the l x A g state occur to l B~2 U i}B^ u ) type 



states via short-axis (long-axis) polarized photons. As per the convention adopted in this 
work, the short (long) axis correspond to y (x) directions. The essential states contributing 
to the linear absorption spectra of various acenes are depicted in [3} The many-particle wave 
functions of the excited states contributing to various peaks in the spectra are presented in 
the EPAPS. 39 While plotting the absorption spectra, we have restricted ourselves to states 
which lie below 6 eV, the estimated value of the ionization potentials of the long acenes ™ 




y-pol. 



Figure 3. Diagram of the essential states involved in the singlet linear optical absorption in 
oligoacenes and their polarization characteristics. The arrows connecting two states imply opti- 
cal absorption, with polarization directions stated next to them. Location of states is not up to 
scale. 



A peak-by-peak detailed description of the computed singlet absorption spectra is pro- 
vided in the supplementary material.^ Here we list the salient featurs of the calculated 
linear optical absorption of octacene, nonacene, and decacene: 

1. Quantitatively speaking, screened parameter spectra are redshifted as compared to 
the standard parameter ones. 

2. Most of the intensity is concentrated in the x-polarized (long- axis polarized) spectra 
originating from the absorption into the 1 B$ U type of states, while the y-polarized 
(short-axis polarized) absorption into the 1 5^ n type states is comparatively weak. 
However, a closer examination reveals that most of the intensity in the x-polarized 
spectrum is derived from the single transition to an n^B^ state (or states which split 
away from it). If we ignore this transition then the short- and long-axis polarized spec- 
tra are of comparable intensity. This aspect of the singlet linear absorption spectrum 
of long acenes is consistent with what is also observed in the shorter acenes.^ 

3. The first peak corresponds to the ^/-polarized transition, to the l 1 ^^ excited state of 
the system. The most important configuration contributing to the many-particle wave 
function of the state corresponds to \H — Y L) excitation, irrespective of the choice of 
the Coulomb parameters employed in the PPP model. 

4. We have observed that the most intense absorption for the oligoacenes is through an x- 
polarized photon to a 1 B 3u state, irrespective of the Coulomb parameters employed in 

10 



the calculations. For acene-n, the many-particle wave function of this state exhibits the 
following general features: (a) for the standard parameter case, single excitations \H — > 
L+n/2 — 1) + c.c, for n = even, and \H — > L + (n — l)/2) +c.c, for n = odd, dominate 
the wave function, while (b) with screened parameters the dominant configurations are 
single excitations \H — > L + n/2) + c.c. for n = even, and \H — > L + {n— l)/2) + c.c. for 
n = odd. The aforesaid difference between the standard and the screened parameters 
is because of different energetic ordering of the symmetries of the molecular orbitals 
for the standard and screened parameters. 



Figure 4. Singlet linear absorption spectra of octacene computed using: (a) standard parameters 
and (b) screened parameters. A uniform line width of 0.1 eV was assumed while plotting the 
spectra. The subscripts attached to the peak labels indicate the polarization directions x and y 
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Figure 5. Singlet linear absorption spectra of nonacene computed using: (a) standard parameters 
and (b) screened parameters. The rest of the information is same as in the caption of E] 
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Figure 6. Singlet linear absorption spectra of decacene computed using: (a) standard parameters 
and (b) screened parameters. The rest of the information is same as in the caption of H] 
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C. Triplet optical absorption calculations 




Figure 7. Diagram of the essential states involved in the triplet absorption spectra of oligoacenes, 
and their polarization characteristics. The arrows connecting two states imply optical absorption, 
with polarization directions stated next to them. Location of states is not up to scale. ISC refers 
to inter-system crossing. 



In optical absorption experiments, one can probe the triplet excited states because fre- 
quently the first singlet excited state Si (l 1 !?^, in the present case) decays to the first triplet 
excited state 7\ (l 3 !?^, in the present case) located below Si, through nonradiative inter- 
system crossing (ISC,), as shown in [71 Once the system reaches the triplet manifold, normal 
optical absorption experiments can be performed to probe higher triplet states. In the 
present work, we restrict ourselves to the triplet one-photon absorption spectra of octacene, 
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nonacene, and decacene from their 1 3 -B^ state, computed using the MRSDCI method. For 
the case of oligoacenes, as per electric-dipole selection rules of the D 2 h point group, the long- 
axis (x-axis) polarized photons cause transitions from the 1 3 -B^ to 3 B± type of states, while 
the short-axis (y-axis) polarized ones lead to the 3 A g type states (cf. Uj. The calculated 
triplet absorption spectra of these oligoacenes are displayed in [8j [9] and]10[ while the wave 
functions of the excited states contributing to various peaks in the spectra are presented 
(cf., Ref. 3 ^). While plotting the triplet absorption spectra, we have been careful to include 
only those states which lie below 6 eV excitation energy (with respect to the l 1 ^" state), 
which is the estimated value of the ionization potentials of the long acenes. 40 A detailed 
description of the characteristics of various peaks in the calculated triplet absorption spec- 
tra of octacene, nonacene, and decacene, is presented in the accompanying supplementary 
materialP^) Below we discuss the salient features of our results: 



1. Similar to the case of singlet absorption, screened parameter spectra are red shifted 
as compared to the standard parameter ones. 



2. Most of the intensity is concentrated in the x-polarized (long- axis polarized) spectra 
corresponding to the aborption into the 3 B± type of states, while the y-polarized 
absorption into the 3 A~ type states is very faint. 



3. From[8j |9j and[TU]it is obvious that the triplet absorption spectrum is dominated by 
two intense x-polarized peaks which are well separated in energy (>2 eV), irrespective 
of the oligoacene in question, or the Coulomb parameters employed in the calculations. 
The first of these peaks is peak I in all the cases, while the second one is either peak 
IV or V, depending upon the oligoacene, or the Coulomb parameters employed. In the 
standard parameter based calculations, peak I is always the second most intense peak, 
while the second of these peaks (IV or V) is the most intense. In the screened parameter 
calculations the situation is exactly the reverse, with peak I being the most intense, 
while peak IV or V being the second most intense peak of the spectrum. Peak I always 
corresponds to the 1 3 B7 excited state of the system, whose wave function is dominated 
by the single excitations \H — ¥ L + 1) + c.c, irrespective of the oligoacene in question, 
or the choice of Coulomb parameters. The second intense peak (IV or V) corresponds 
to a higher 3 B± type state of acene-n, whose wave function mainly consists of: (a) 
for the standard parameters, double excitations \H — >■ L; H — (n/2 — 1) — > L) + c.c. 
for n = even, and \H — >■ L; H — (n — l)/2 — ¥ L) + c.c. for the n = odd case, while 
(b) for the screened parameters, double excitations \H — >■ L; H — n/2 — > L) + c.c. 
for n = even, and \H — >■ L; H — (n — l)/2 — > L) + c.c. for the n = odd case. This 
difference between the standard and the screened parameter results is because of the 
different energetic ordering of the molecular orbital symmetries in longer acenes, for 
the two sets of Coulomb parameters. 



4. Other peaks in the spectrum correspond to either x or y polarized transitions to the 
higher excited states of the system, which are described in detail (cf., Ref™' 1 
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Figure 8. Triplet absorption spectra of octacene from the 1 3 -B^ state computed using: (a) standard 
parameters and, (b) screened parameters. A uniform line width of 0.1 eV was assumed while plotting 
the spectra. The subscripts attached to the peak labels indicate the polarization directions x and 
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Figure 9. Triplet absorption spectra of nonacene from the l^B^ state computed using: (a) standard 
parameters and, (b) screened parameters. The rest of the information is same as in the caption of 
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Figure 10. Triplet absorption spectra of decacene from the 1 B^ u state computed using: (a) the 
standard parameters and, (b) the screened parameters. The rest of the information is same as in 
the caption of [8] 
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D. Comparison of singlet and triplet absorption calculations of octacene and 
nonacene with the experimental results 

In [V] the results of the experimental absorption of octacene and nonacene, reported by 
Tonshoff and Bettinger,^ have been presented, showing the excitation energies, relative in- 
tensities, and tentative assignments of the UV/vis electronic absorptions. Next, we compare 
our computed singlet and triplet absorption results, with the expermental ones, in terms of 
the peak energies (in eV) and relative intensities/oscillator strength. The relative oscillator 
strength (ROS) in the computed absorption spectra is the ratio of the calculated oscillator 
strength of a given state, with respect to that of the most intense state. The relative os- 
cillator strengths of the peaks in the calculated absorption spectra are compared with the 
reported relative intensities (RI) of the corresponding peaks in the experiment spectra. 

1. Octacene 

Comparing the experimental values of the energies of the absorption peak of octacene 
from the [Va, with our calculated values for the singlet case (cf., RefP*'), the lowest state 
Si (1.54 - 1.86 eV, RI= 0.025 for the 1.54 eV) of the experimental results matches nicely with 
the first peak at 1.49 eV (ROS= 0.050) obtained in the screened parameter calculations, than 
the corresponding one for the standard parameter case located at 2.24 eV. The second state 
of the experimental results S2 (2.40 eV, RI=0.02) again agrees much better with the screened 
parameter result, 2.65 eV (ROS=0.047) than with the standard parameter one located at 
3.34 eV. The third state of the experimental results, S3 (2.54 eV, RI=0.050) has a better 
agreement with the third state of the screened parameter spectrum 2.97 eV (ROS=0.046), 
than with corresponding standard parameter peak at 4.17 eV. The most intense peak in the 
experimental results corresponds to S-j (3.68-3.78 eV, RI=1.0) lies in between the highest 
peak of the screened and standard parameter calculations located at 3.38 eV (ROS=1.0), 
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and 4.17 eV (ROS=1.0), respectively. We note that with the increasing peak energy, the 
agreement between the experiments and theory begins to deteriorate. Nevertheless, even in 
the worst case scenario of the most intense peak (£7 of the experiments), the disagreement 
between the screened parameter results and the experiments is below 10%, as far as the peak 
location is concerned. This disagreement, which is fairly acceptable from a quantitative 
aspect, is possibly due to the reduced computational accuracy of our approach because 
of the high energies of the states involved. As far as the comparison of the experimental 
results with the triplet absorption calculations (cf., RefJ^j) is concerned, the first peaks in 
the computed triplet spectra, 3.04 eV (ROS=0.839) for the standard parameter, and 1.97 
eV (ROS=1.0) for the screened parameters, disagree completely with the location of Si, as 
well as its relative intensity in the experimental spectrum. The same trend holds for the 
higher excited states as well. Thus, we conclude that the experimental absorption spectrum 
of octacene is indeed from the 1 X A~ ground state of the system, confirming yet again that 
the ground state of octacene has the singlet multiplicity. 

We also note that Raghu et alW^ predicted a very large optical gap of 2.60 eV for octacene, 
based upon their DMRG calculations, employing the PPP Hamiltonian. 



Table V. The experimental values of the wavelenj 
of the UV/vis electronic absorptions of octacene 



;ths, relative intensities and tentative assignments 
(left) and nonacene (right j^ 



E(eV) 


I (rel.) 


Octacene State 








1.54 


0.025 


Si 


E(eV) 


I (rel.) 


Nonacene State 


1.69 


0.011 


1.43 


0.020 


£1 


1.73 


0.016 


1.58 


0.006 


1.86 


0.007 


1.62 


0.016 


2.40 


0.02 


s 2 


2.33 


0.033 


S3 


2.54 


0.05 


£3 


2.50 


0.023 


£4 


2.72 


0.05 


2.67 


0.030 


£5 


3.16 


0.08 


£4 


2.80 


0.048 


Se 


3.29 


0.49 


£ 5 


2.97 


0.19 


s 7 


3.33 


0.43 


3.14 


0.52 


3.48 


0.46 


£ 6 


3.42 


0.14 


Ss 


3.68 


0.90 


£7 


3.66 


1.00 


£9 


3.78 


1.00 









2. Nonacene 



Similarly, comparing the experimental values of the energies of the absorption peak of 
nonacene from the|Vp, with our calculated values for the singlet case (cf., RefP*'), the lowest 
state £1 (1.43 - 1.62 eV, RI =0.020 for 1.43 eV) of the experimental results matches nicely 
with the first peak of 1.46 eV (ROS = 0.051) obtained in the screened parameter calculations, 
than the corresponding one for the standard parameter case located at 1.82 eV. The second 
state of the experimental results, £ 3 (2.33 eV, RI = 0.033) again agrees much better with the 
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second peak of the screened parameter result, 2.45 eV (ROS = 0.043) than with the standard 
parameter one located at 2.79 eV. The third state of the experimental results, £4 (2.50 eV, 
RI = 0.023) has a better agreement with the third state of the screened parameter spectrum 
2.71 eV (ROS = 0.020) than with the corresponding standard parameter peak at 3.80 eV. 
The most intense peak in the experimental results is Sg (3.66 eV, RI = 1.0) lies in between 
the highest peak of the screened and standard parameter calculations located at 3.32 eV 
(ROS = 1.0), and 3.80 eV (ROS = 1.0), respectively. For this most intense peak, the stan- 
dard parameter results appear to agree slightly better with the experiments, as compared 
to the screened ones. But, keeping in mind the lower energy peaks discussed above, overall 
the screened parameter based results have a much better agreement with the experiments, 
just as in case of octacene. Similar to octacene, we again note that with the increasing peak 
energy, the agreement between the experiments and theory begins to deteriorate, which we 
again attribute to the reduced computational accuracy for higher energies. As far as the 
comparison of the experimental results with the triplet absorption calculations (cf., RefP^) 
is concerned, the first peaks in the computed triplet spectra, 2.56 eV (ROS=0.785) for the 
standard parameter, and 1.86 eV (ROS=1.0) for the screened parameters, disagree com- 
pletely with the location of Si, as well as its relative intensity in the experimental spectrum. 
The same trend holds for the higher excited states as well. Thus, we conclude that the 
experimental absorption spectrum of nonacene is indeed from the ^- l A~ ground state of the 
system, confirming yet again that the ground state of nonacene has the singlet multiplicity. 
Again, we note that the DMRG based calculations employing the PPP model, performed 
by Raghu et al^ predict an unrealistically large optical gap of 2.59 eV for nonacene. Our 
values of the calculated optical gaps of octacene, nonacene, and decacene are in very good 
agreement with the estimated optical gap of polyacene 1.18 ±0.06 eV, reported by Tonshoff 
and BettingerP^ 

To conclude, the experimentalists are confident that the absorptions which they observed 
are all from the ground states of octacene and nonacene. Because, our singlet absorption 
spectra (as against the triplet absorption) computed using the screened parameters are in 
excellent agreement with those experiments, we conclude that the ground states of octacene 
and nonacene are of singlet multiplicity (l 1 ^), in perfect agreement with our singlet-triplet 
crossover calculations. Although the experiments have not been performed on acenes longer 
than nonacene, the trend emerging from our calculations leads us to conclude that the 
ground state will be of singlet type in those systems as well. 



IV. CONCLUSIONS 

To summarize, we presented large-scale MRSDCI calculations on the electronic structure 
and optical properties of oligoacenes, with focus on the longer acenes, namely, octacene, 
nonacene, and decacene. By performing such calculations on the lowest singlet and triplet 
states of oligomers ranging from naphthalene up to decacene, we established that the ground 
state in oligoacenes has singlet multiplicity, with a singlet-triplet gap of approximately 0.5 eV 
even for decacene. The trends visible from our calculations rule out a singlet-triplet crossover 
for the ground states of longer oligoacenes as well. This result of ours has thus resolved an 
old speculation predicting that nonacene onwards, the ground state of oligoacenes will be of 
triplet multiplicity.^ 

Moreover, the many-body wavefunction analysis of the correlated singlet ground state 
l 1 ^!" reveals increasing contribution of configurations with two open shells, accompanied 
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with the decreasing one from the closed-shell Hartree-Fock reference state, with the increas- 
ing chain length. Thus our calculations predict an open-shell diradical character for the 
singlet ground state of longer acenesP^ 

As far as the singlet linear optical absorption is concerned, in all the acenes, the first peak 
is due to a y-polarized transition to the l 1 -^^ state, corresponding to the HOMO to LUMO 
transition. The most intense state is the x-polarized transition to a 1 B^ U state, which is also 
dominated by single excitations. When we compare our singlet linear absorption spectra 
of octacene and nonacene with the experimental ones ™ excellent agreement is obtained on 
the important peak locations and intensity profiles. Furthermore, the measured absorption 
spectra of longer acenes show no resemblance with our computed triplet absorption spectra, 
confirming once again the conclusion that the ground state of the longer acenes is indeed 
singlet in nature. 

Our calculations on the one-photon triplet absorption spectra predict two intense x- 
polarized absorptions, which are well separated in energy. Besides these two, there are a 
number of weaker peaks which are either x or y polarized. This is in sharp contrast to the 
singlet absorption which predicts only one intense peak. The existence of two well separated 
x-polarized peaks in the triplet absorption spectrum is one of the most important predictions 
of this work, and can be tested in future experiments on oriented samples of longer acenes. 

We also performed singlet and triplet optical absorption calculations on decacene, a 
molecule which has not been synthesized yet. We are hopeful that in future, once decacene 
is synthesized in the laboratory, our theoretical predictions could be tested in experiments. 

In this paper we restricted ourselves to the low-lying excited states of longer acenes 
which contribute to their linear optical properties. However, not many calculations have 
been performed as far as the nonlinear optical properties of these materials are concerned. 
In particular, it will be of interest to compute the nonlinear susceptibilities corresponding 
to two-photon absorption, and third harmonic generation. Both these nonlinear optical 
processes have the capability to probe the higher excited states of polyacenes, which is 
essential in order to obtain a deeper understanding of the optical response of it electrons. 
At present, studies along these directions are underway in our group. 
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Supplementary Information: Detailed 
Description of Calculated Spectra 



V. SINGLET LINEAR ABSORPTION 

Here we present a detailed description of the calculated singlet linear absorption spectra 
of octacene, nonacene, and decacene, presented in Figures 4, 5 and 6 of the main text. Our 
calculations reveal the following features: 

1. For all the oligoacene, the first peak is due to a y-polarized transition, to the l 1 -^^ 
excited state of the system, whose wave function is dominated by the \H — >■ L) single 
excitation, irrespective of the choice of the Coulomb parameters employed in the PPP 
model. 

2. The second peak also corresponds to a y-polarized transition, to the 2 1 5^ i excited state 
of the system. The most important configuration contributing to the many-particle 
wave function of this state is \H — 1 — > L + l)excitation, irrespective of the choice of 
the Coulomb parameters. 

3. The nature of the third peak is dependent upon the Coulomb parameters employed in 
the PPP model. For the standard parameter case, this peak always corresponds to the 
re-polarized, ^B^ u excited state, signalling the onset of the most intense absorption 
feature in the system. For octacene, it is the single most intense peak, whose wave 
function mainly consists of single excitations \H — > L + 3) + c.c. (c.c. denotes the 
charge conjugated configuration). For nonacene and decacene, perhaps due to band 
formation, it is the first of the two adjacent intense peaks which are x polarized. For 
nonacene it is the most intense peak, and the corresponding many-body wave function 
is dominated by the single excitation \H — > L + 4) + c.c. For decacene, however, it is 
the second most intense feature, and the double excitations \H — >■ L; H — Y L + l) + c.c. 
contribute the most to its many-body wave function. 

In the screened parameter calculations, the third peak is a faint peak, containing 
a mixture of x and y polarized transition to the states, l 1 ^^ and ^B^ of the system. 
For all the oligoacenes, the double excitation \H — > L; H — > L + 1} + c.c. contributes 
the most to the l 1 -^^ state, and the single excitation \H — >■ L + 2) + c.c. dominates 
the ^B^u wave function. 

4. For the standard parameter case, the fourth peak is a faint peak for octacene, con- 
sisting of a mixture of x- and ^/-polarized transitions to the states, 2 1 B^ U and 4 1 5^ u . 
Double excitations \H — > L; H — > L + 1} +c.c. contribute the most to the 2 1 B^ U state, 
and single excitations \H — 2 — > L + 2) + c.c. to the 4 1 5^ i state. For nonacene and de- 
cacene, however, it is an intense re-polarized feature corresponding to the state 2 1 B^ U , 
which is adjacent to their l^-B-^ state mentioned above. For the case of decacene, it 
also the most intense absorption of the system. Both for nonacene and decacene, the 
single excitations \H — > L + 4) + c.c. contribute the most to the many-particle wave 
function of this state. 



20 



For the screened parameter calculations, the fourth peak corresponds to the most 
intense absorption, through an a;-polarized photon, to the 2 l B^ u state, in case of 
octacene and nonacene. For decacene, however, it appears as a shoulder to the most 
intense peak, consisting of a mixture of x and y polarized transitions, to the 2 X B^ U and 
4 1 Z?^ U states. The many-particle wave function of the 2 X B^ U state for octacene and 
nonacene is dominated by the single excitations \H — > L + 4) + cc, while for decacene 
the double excitations \H — >■ L + 1; H — > L + 2) + c.c. dominate. The 4 1 5^ u state of 
decacene has the maximum contribution from the single excitation \H — 2 — >■ L + 2). 

5. For the standard parameter case, the fifth peak is a very faint feature consisting of 
y-polarized transition to the state l^B^ for octacene, whose wave function is domi- 
nated by the single excitation \H — 2 — > L + 2). However, for nonacene and decacene, 
it corresponds to the mixture of x and y polarized transitions to the states, 3 l B^ u and 
^B^u, respectively. The double excitations \H — > L + 1; H — 1— > L + l) + c.c. for the 
^B^ state, and the triple excitation \H — > L; H — > L; H — 1 — > L + 1} for the ^B^ 
state, contribute the most to their wave functions. 

For the screened parameter case, it is a faint peak consisting of a mixture of x and y 
polarized transitions to the states, ^B^ and 4 1 5^ t , respectively, for both octacene 
and nonacene. For decacene, however, it is the most intense peak corresponding to x 
polarized transition to the state 4 1 _B 3 ' U . The most important configuration contribut- 
ing to the many-particle wave function of the state for octacene and nonacene, is a 
mixture of double excitations \H — >■ L + l; H — > L + 2) +c.c. for 3 l B^ u state, and single 
excitations \H — 2 — > L + 2) for the 4 1 5j t state. For decacene, the single excitations 
\H — > L + 5) + c.c. contribute the most to the wavefunction of this state. 

6. For the standard parameter case, the sixth peak corresponds to a mixture of x and 
y polarized transitions to the states, 4 1 Z? 3 ' M and ^B^, respectively, for octacene. 
The many-body wave functions of both these states derive predominant contribu- 
tions from double excitations: \H — > L; H — 1 — > L + 2) + c.c. for the 4 1 _B 3 I U , and 
\H — >■ L + 1; H — ¥ L + 3) + c.c. for the ^B^ state. For nonacene (decacene), it 
corresponds to a ^/-polarized transition to the state lO 1 !?^ (9 1 5^ u ), whose many-body 
wave function derives maximum contribution from the single excitation \H— 2 — > L+2). 

The same peak with the screened parameters is a mixture of x and y polarized transi- 
tions for the case of octacene and nonacene. For octacene, transitions to states Q l B^ u 
and 8 l B2 U constitute this peak, with double excitations \H — > L; H — > L + 1) + c.c. 
and \H — > L; H — > L + 5) + c.c, respectively, contributing most to their wave func- 
tions. For nonacene, this peak consists of three states 4 1 5^ u and Q 1 B2 U and l^B^ 
with double excitation \H — > L; H — 1 — > L + 2) + c.c. contributing the most to the 
4 1 i?^ u , triple excitation \H — >■ L; H — >■ L; H — 1 — > L + 1) to the 6 1 B 2 \ t , and the single 
excitation \H — 2 — > L + 2) to the ^B^ state. However, for decacene, the peak is 
purely x polarized, due to transition to the state Q 1 B^ U whose wave function derives 
maximum contribution from the double excitations \H — > L; H — y L + 1) + c.c. 

7. The seventh peak in the standard parameter spectrum is both x and y polarized for 
all the three oligomers. For octacene, the peak involves 5 1 B^ U and 10 1 B 2 + U states, with 
single excitations \H — > L + 7) + c.c. and \H — 2 — ¥ L + 2), respectively, contributing 
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the most to their wave functions. For nonacene states 7 1 B% U and 12 1 5^ u constitute the 
peak, with single excitation \H — > L + 8)+c.c. providing the maximum contribution 
to the 7 1 B^ U state, and triple excitation \H — > L; H — 1 — > L + 1; H — 1 — > L + 1) 
to the 12 1 i?^ u state. For decacene, the peak involves Q X B^ U and ll 1 !?^ states, 
with double excitations \H — > L; H — 1 — > L + 2) + c.c. and triple excitations 
\H — > L; H — 1 — >■ L + 1; H — 1 — >■ L + 1), respectively, contributing the most to their 
wave functions. 

With the screened parameters, the seventh peak of octacene and nonacene is x polar- 
ized, while for decacene it has mixed x and y polarizations. Both for octacene and 
nonacene states 8 1 B^ U constitute this peak, with single excitation \H — 1 —¥ L + 5) + c.c. 
and \H — 1 — > L + 6} + c.c. providing main contributions to the wave functions of oc- 
tacene, and nonacene, respectively. However, for decacene, transitions to states 8 1 B^ U 
and lO 1 !?^ form the peak, with double excitations \H — > L; H — 1 — > L + 4) + c.c. 
and single excitations \H — 1 — > L + 7) + c.c, respectively, contributing the most to 
their wave functions. 

8. With the standard parameters, the eighth peak does not exist for octacene, while 
for nonacene and decacene, it has mixed x and y polarizations. For nonacene, the 
states constituting the peak are 8 1 B^ U and lG 1 !^, with single excitations \H — 1 — >■ 
L + 5) + c.c. and triple excitations \H — > L; H — > L + 2; H — 1 — > L + 1) + c.c, 
respectively, contributing to their wave functions. For decacene, three states 7 1 B% U , 
1A 1 B2 U , and lS 1 !?^, form this peak, with single excitations \H — > L + 8) + c.c, 
\H — > L + 12) + c.c, and \H — > L + 9} + c.c, respectively, contributing the most to 
their many-particle wave functions. 

With the screened parameters case, the eighth peak is x and y polarized for oc- 
tacene and decacene, while it is only x polarized for nonacene. For octacene, states 
Q 1 B^ U and l&B^ form this peak, whose wave functions, respectively, derive maximum 
contributions from single excitations \H — 2 — > L + 4) + c.c, and triple excitations 
\H — > L; H — 1 — > L + 1; H — 1 — > L + 1). In case of nonacene, the transition to 
the state Q l B^ u leads to this peak, with single excitations \H — 2 — > L + 4) + c.c. 
contributing the most to its many-particle wave function. For decacene, states Q 1 B^ U 
and 14 l B2 U constitute the peak, with the single excitations \H — > L + 8) + c.c. and 
\H — > L + 9) + c.c, respectively, providing the largest contributions to their wave 
functions. 

9. The ninth peak, for the standard parameter case, exists only for decacene, and has 
a mixed x and y polarization, with states 8 l B% u and 17 1 B2 U forming the peak. The 
single excitations \H — 2 — > L + 4) + c.c. contribute the most to S 1 ^^ state, while 
double excitations \H — > L; H — 1 — > L + 4) + c.c. to the 17 l B2 U state. 

With the screened parameters, the ninth peak does not exist for nonacene, and it is 
x and y polarized for octacene and decacene. For octacene, states 12 l B^ u and l&B^u 
constitute this peak, with the double excitations \H — > L; H — 1 — > L + 6) + c.c. and 
single excitations \H — 1 — >■ L + 7) + c.c, respectively, contributing the most to their 
wave functions. For decacene, the states ll 1 !?^ and 16 1 B2 U form this peak, with the 
single excitations \H — 2 — )■ L + 5) + c.c. and \H — 1 — )■ L + 7) + c.c, respectively, 
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providing the maximum contributions to their wave functions. 

Table VI: Excited states contributing to the singlet linear ab- 
sorption spectrum of octacene computed using the MRSDCI 
method coupled with the standard parameters in the PPP 
model Hamiltonian. The table includes many particle wave 
functions, excitation energies, dipole matrix elements, and 
relative oscillator strengths (ROS) of various states with re- 
spect to the l l A g ground state. DF corresponds to dipole 
forbidden state. Below, i +c.c.'' indicates that the coefficient 
of charge conjugate of a given configuration has the same sign, 
while '—c.c.' implies that the two coefficients have opposite 





signs. 








Peak 


State 


E (eV) 


Transition 
Dipole (A) 


ROS 


Wave Functions 


DF 


l 1 ^ 


2.31 








\H -+ L; H -+ L + 1 ) - c.c.(0.4948) 
\H -> L;H -1-+ L + 2) - c.c.(0.1938) 


I 


1^ 2 + 


2.24 


0.905 


0.034 


\H->L) (-0.8471) 
\H -»■ L;H -»• L; # - 1 -> L + l) (0.0948) 


II 


*BL 


3.34 


0.641 


0.025 


\H-1-+L+1) (-0.5826) 
\H -> L + 2) - c.c. (0.4284) 


III 


1^3+ 


4.17 


3.622 


1.000 


\H ^ L + 3) +c.c.(0.4622) 
\H -+ L;H -> L + l) + c.c.(0.2614) 


IV 


4^+ 


4.51 


0.367 


0.011 


\H -2-> L + 2) (-0.4148) 
\H-l-+L + 4) + c.c.(0.4020) 




2^3+ 


4.57 


1.079 


0.097 


\H -+ L;H -+ L + l) + c.c.(0.3214) 
\H ^ L + 3) + c.c. (0.3214) 


V 


7'Bl 


5.46 


0.375 


0.014 


\H -2-> L + 2) (0.3574) 
\H -»■ L + 6) - c.c. (0.3489) 


VI 


9^ 2 + 


5.64 


0.261 


0.007 


\H-+L + l;H-+L + 3) + c.c.(0.3883) 
\H -+ L;H -+ L + 5) + c.c.(0.2768) 




4^3+ 


5.65 


0.335 


0.012 


\H -+ L;H -1-+ L + 2) - c.c.(0.4696) 
\H -+ L; H -+ L + 1 ) + c.c.(0.1289) 


VII 


10^2+ 


5.89 


0.203 


0.004 


\H -2-+ L + 2) (-0.3657) 
\H -3-» L + 3) (0.3074) 




5^3+ 


5.94 


0.307 


0.010 


\H -+ L + 7) - c.c. (0.3397) 

\H -+ L;H -+ L + 4) + c.c.(0.2494) 



Table VII: Excited states contributing to the singlet linear ab- 
sorption spectrum of octacene computed using the MRSDCI 
method coupled with the screened parameters in the PPP 
model Hamiltonian. The rest of the information is same as 
that in NT 



23 



Peak 


State 


E(eV) 


Transition ROS 
Dipole (A) 


Wave Functions 


DF 


l 1 ^ 


1.59 





\H -»■ L;H -+ L + l) - c.c.(0.5078) 
\H -»• L;H - 1 ->■ L + 2 ) - c.c.(0.1949) 


I 


1'Bl 


1.49 


1.241 0.050 


\H-*L) (0.8503) 
\H -+ L;H -+ L;H - 1 -> L + l) (0.0929) 


II 


I'Bl 


2.65 


0.897 0.047 


\H-1-+L + 1) (-0.7244) 
\H -+ L + 2) + c.c.(0.2638) 


III 


Z'Bl 
1^3+ 


2.89 
2.97 


0.440 0.012 
0.845 0.046 


|#^L + 2) +c.c.(0.5326) 

\H-1->L + 1) (-0.3290) 

\H-*L]H->L + 1) - c.c.(0.4942) 

\H->L+1;H-1->L+1) - c.c.(0.2522) 


IV 


I'Bl 


3.38 


3.675 1.000 


\H -+ L + 4) +c.c.(0.5831) 
|LT-l->-L + 5) + c.c. (0.1004) 


V 


Z'Bl 
^Bl 


3.91 
3.97 


0.410 0.014 
0.641 0.036 


|#^L + l;F^L + 2) - c.c.(0.4028) 

jiT-).L + l;iT->.L + 4) - c.c.(0.3003) 

|#-2^L + 2) (-0.5280) 

|LT - 1 -+ £ + 3 ) + c.c.(0.3722) 


VI 




4.52 

4.57 


0.230 0.005 
0.403 0.016 


\H ^ L;H ^ L + 5) - c.c.(0.4706) 

|i?^L;#-l^L + 4) - c.c.(0.3075) 

|JET -»■ L;H -*• L + 1 ) - c.c.(0.3419) 

\H ^ L;H ^ L + 3) - c.c.(0.2890) 


VII 


8^3 + . 


4.80 


1.301 0.178 


|iT - 1 -+ L + 5 ) + c.c.(0.5537) 
\H -+ L;H - 1 -> L + 2) - c.c.(0.1321) 


VIII 


9^3 + . 

13^+ 


5.14 
5.22 


0.410 0.019 
0.418 0.020 


|# - 2 ^ L + 4 ) + c.c.(0.5720) 

|Lf^L + 2;F^L + 3) - c.c.(0.0943) 

\H-+L;H-1-*L + 1;H-1-*L + 1) (0.3505) 

\H - 1 -> L + 7 ) + c.c.(0.2709) 


IX 




5.53 

5.66 


0.440 0.023 
0.154 0.003 


|LT - 1 -»• L + 7 ) + c.c.(0.3983) 

|Lf-3^L + 3) (0.3636) 

\H -+ L; Lf - 1 -+ L + 6) - c.c.(0.4014) 

|L r ^L + 2;F-2^L + l)- c.c.(0.3019) 



Table VIII: Excited states contributing to the singlet linear 
absorption spectrum of nonacene computed using the MRS- 
DCI method coupled with the standard parameters in the 
PPP model Hamiltonian. The rest of the information is same 







as tha 


t in 


VI 






Peak 


State 


E (eV) 


Transition ROS 
Dipole (A) 


Wave Functions 


DF 


l 1 ^ 


1.86 





\H ^ L;H ^ L + l) + c.c.(0.4872) 
\H -> L + 1;H -> L + 2) + c.c.(0.2027) 


I 


1^ 2 + 


1.82 


1.328 0.092 


\H->L) (0.8290) 
\H -> L;H -> L;H -l-> L + l) (0.1066) 
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II 


2^ 2 + 


2.79 


0.733 0.043 


\H-l-+L + l) (0.5506) 
\H -+ L + 2) +c.c.(0.4342) 


III 


1^3+ 


3.80 


3.037 1.000 


|tf ^L + 4) - c.c. (0.3757) 

|if ->■ L; H ->■ L + 1 ) - c.c.(0.3274) 


IV 


2^3+ 


4.12 


2.583 0.786 


|#^L + 4) - c.c. (0.3943) 

\H ^ L;H ^ L + l) - c.c.(0.2791) 


V 


S 1 ^ 


4.62 
4.63 


0.439 0.025 
0.549 0.040 


|#^L;#^L;#-1^L + 1) (0.6215) 

\H-1-+L + 1) (0.2357) 

|if^-L+l;iT-l^L+l)- c.c.(0.3339) 

|fl'-)-L + l;fl'-)-L + 2) - c.c.(0.3035) 


VI 


10^2+ 


5.29 


0.457 0.032 


\H -2-> L + 2) (-0.3824) 
\H->L + l;H->L + 4,) -fee. (0.2681) 


VII 


7^3+ 


5.71 
5.72 


0.277 0.013 
0.318 0.017 


\H -+ L + 8) - c.c. (0.3021) 

\H -+ L;H -+ L + 3) + c.c.(0.2619) 
|il ->• -L; H - 1 ->• L + 1; H - 1 ->• L + 1 ) (-0.3293) 

\H -> L;H -l-> L + 4) - c.c.(0.2869) 


VIII 


16^2+ 
8^3+ 


5.94 
5.95 


0.287 0.014 
0.996 0.169 


\H -> L; H -> L + 2; H - 1 -> L + 1) - c.c.(0.3742) 

\H -+ L; ff - 1 ->• L + 4) - c.c.(0.2500) 

|# - 1 -»• £ + 5 ) - c.c.(0.3932) 

\H -2^ L + 4) - c.c.(0.2993) 



Table IX: Excited states contributing to the singlet linear ab- 
sorption spectrum of nonacene computed using the MRSDCI 
method coupled with the screened parameters in the PPP 
model Hamiltonian. The rest of the information is same as 





that in 


V 


[ 




Peak 


State 


E (eV) 


Transition ROS 
Dipole (A) 


Wave Functions 


DF 


l 1 ^ 


1.51 





\H -+ L; H -» L + 1 ) + c.c.(0.5143) 

\H -> L;H -1-+ L + 2) - c.c.(0.2089) 


I 


1^2+ 


1.46 


1.316 0.051 


\H -+L) (0.8551) 
\H -+ L; H -+ L; H - 1 -+ L + 1 ) (0.1017) 


II 


2^ 


2.45 


0.935 0.043 


\H-l-±L + l) (0.7260) 
|F^L + 2) -c.c.(0.2689) 


III 


l 1 ^ 
3^2+ 


2.71 

2.77 


0.611 0.020 
0.507 0.014 


\H -»• L; # -»■ L + 1 ) - c.c.(0.4901) 

|F^L + l;i?-l^L + l) - c.c.(0.2454) 

|F^L + 2) -c.c.(0.5434) 

|iI-l-)-L + l) (-0.3462) 


IV 


2^ 


3.32 


3.887 1.000 


|F^L + 4) -c.c.(0.5689) 
\H-l^L + 6) +c.c.(0.1185) 


V 


3^3 + « 

4^2+ 


3.63 
3.69 


0.782 0.044 
0.559 0.023 


\H->L+l;H->L + 2) + c.c.(0.3441) 

\H-1->L + 1;H-+L + 1) - c.c.(0.3221) 

\H -2-> L + 2) (-0.4715) 

\H-l-+L + 3) + c.c.(0.4395) 
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VI 


4^3+ 


3.93 


0.369 


0.011 


\H -> L;H - 1 -» L + 2) - c.c.(0.4953) 
\H -+ L;H -+ L + 3) - c.c.(0.2163) 




e 1 ^ 


4.01 


0.690 


0.038 


\H -+ L; H -+ L; il - 1 ->• L + 1 ) (0.5954) 
|F^L + 5) +c.c.(0.3166) 




7'Bl u 


4.12 


0.692 


0.039 


\H -2^ L + 2) (-0.5949) 
|#-l^L + 3) + c.c.(0.3994) 


VII 


#B& 


4.61 


1.450 


0.193 


|#-l^L + 6) + c.c.(0.5262) 
\H -+ L;H - 1 ->• L + 2) - c.c.(0.1473) 


VIII 


^Bl 


4.97 


0.338 


0.011 


|#-2^L + 4) + c.c.(0.5823) 
|F^L + 8) +c.c.(0.0699) 



Table X: Excited states contributing to the singlet linear ab- 
sorption spectrum of decacene computed using the MRSDCI 
method coupled with the standard parameters in the PPP 
model Hamiltonian. The rest of the information is same as 







that i] 


i VI 






Peak 


State 


E (eV) 


Transition ROS 
Dipole (A) 


Wave Functions 


DF 


l 1 ^. 


1.72 





\H -» L;H -+ L + 1 ) - c.c.(0.4790) 
\H->L + l;H->L + 2) + c.c.(0.2119) 


I 


1^ 2 + 


1.79 


1.423 0.084 


\H->L) (-0.8203) 
\H-1-+L + 1) (-0.1213) 


II 


2^+ 


2.64 


0.776 0.037 


\H-1->L+1) (-0.5607) 
\H -> L + 2) - c.c. (0.4239) 


III 


1^3+ 


3.68 


2.554 0.553 


\H -+ L;H -+ L + l) + c.c.(0.3710) 
\H -> L + A) +c.c.(0.2865) 


IV 


2^3+ 


4.00 


3.293 1.000 


\H -+ L + A) + c.c. (0.4451) 
\H-+L + l]H-+L + 2) - c.c.(0.2322) 


V 


5^ 2 + 
3^+ 


4.40 
4.46 


0.505 0.026 
0.739 0.056 


\H -+ L; H -> L; H - 1 -+ L + 1 ) (-0.6203) 

\H-l-+L + l) (-0.2380) 

\H ->L + 1;H-1-+L + 1) + c.c.(0.3459) 

\H->L + l;H-*L + 2) - c.c.(0.2800) 


VI 


9^ 2 + 


5.05 


0.392 0.018 


\H -2-+ L + 2) (-0.3901) 
\H-l-+L + 3) (-0.2482) 


VII 


6^3+ 


5.45 
5.45 


0.337 0.014 
0.351 0.015 


\H -> L;H - 1 -» L + 2) - c.c.(0.3992) 

\H ^ L + 8) - c.c. (0.2341) 

\H-+L;H-1-+L + 1;H-1-+L + 1) (-0.4912) 

\H -+ L; H -»• L; iJ - 1 ->• L + 1 ) (-0.2510) 


VIII 


^Bt 
14^+ 

^Btu 


5.59 
5.62 
5.63 


0.584 0.044 
0.332 0.014 
0.326 0.014 


|#^L + 8) (0.3796) 
\H-l-+L + 6) - c.c.(0.2866) 

\H ^ L + 12) +c.c.(0.4160) 
\H - 1 -+ L + 7 ) + c.c.(0.2309) 

\H -+ L + 9) +c.c.(0.3036) 
\H-l->L + 7)+ c.c.(0.2945) 
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IX 


8^3+ 


5.89 


0.913 


0.113 


\H -2^ L + 4) - c.c.(0.3596) 
\H-l-+L + 6) - c.c.(0.2949) 




17^+ 


5.89 


0.149 


0.003 


\H->L;H-l->L + 4) - c.c.(0.4687) 
\H -> L; H -»• L; H - 2 -»• L + 2) (-0.1847) 



Table XI: Excited states contributing to the singlet linear ab- 
sorption spectrum of decacene computed using the MRSDCI 
method coupled with the screened parameters in the PPP 
model Hamiltonian. The rest of the information is same as 





that in 


VI 








Peak 


State 


E (eV) 


Transition 
Dipole (A) 


ROS 


Wave Functions 


DF 


l 1 ^ 


1.15 








\H -» L;H ->• L + 1 ) - c.c.(0.4851) 
\H->L+l]H-l-*L + 2) - c.c.(0.2102) 


I 


I 1 **, 


1.27 


1.369 


0.067 


\H-*L) (-0.8338) 
\H -> L; H -»• L; H - 1 -»• L + 1 ) (+0.1112) 


II 


2^ 2 + 


2.16 


0.929 


0.052 


|iJ-l->L+l) (+0.7110) 
#^L + 2) + c.c.(0.2625) 


III 


3^ 2 + 


2.42 


0.417 


0.012 


|#^L + 2) +c.c.(0.5220) 
iJ-l^-L+1) (-0.3185) 




^ 


2.42 


0.295 


0.006 


\H -» L; F -»• L + l) + c.c.(0.4762) 
|iJ-l-^L;i!-l-)-L + l)+ c.c.(0.2590) 


IV 


2^3 + U 


3.20 


1.872 


0.315 


|iI-^L + l;iT^L + 2)+ c.c.(0.3455) 
|#^L + 5) +c.c.(0.2733) 




4^ 


3.21 


0.706 


0.045 


|#-2^L + 2) (0.5140) 
\H - 1 -»• L + 3 ) + c.c.(0.3975) 


V 


4^ 


3.35 


3.259 


1.000 


|#^L + 5) +c.c.(0.4970) 
\H ^ L;H -1^ L + 2) - c.c.(0.2002) 


VI 


e 1 ^ 


3.81 


0.273 


0.008 


\H -*■ L;H -»• L + 1 ) - c.c.(0.3703) 
\H ^ L;H ^ L + 3) + c.c.(0.3689) 


VII 


s 1 ^ 


4.32 


0.184 


0.004 


\H ^ L;H - 1 -^L + 4) -c.c.(0.2891) 

|#^L + 6;F-5^L + 2)- c.c.(0.2750) 




io 1 ^ 


4.33 


0.478 


0.028 


|iJ- 1 -^L + 7) -c.c.(0.2951) 
iJ - 2 ->■ L + 4) - c.c.(0.2793) 


VIII 


S 1 ^ 


4.51 


0.330 


0.014 


|#^L + 8) +c.c.(0.4467) 
\H -> L + 2;H - 1 ->■ L + 2) + c.c.(0.2155) 




14^+ 


4.56 


0.376 


0.018 


|#^L + 9) +c.c.(0.4944) 
#-3^L + 7) (-0.2478) 


IX 


11^3+ 


4.77 


0.754 


0.076 


|if-2^L + 5)+ c.c.(0.5543) 
|#^L + 3;F-1^L + 1)- c.c.(0.2076) 




^bL 


4.77 


0.520 


0.036 


|iT- 1 -»• £ + 7) - c.c.(0.3795) 
|#-3^L + 7) (0.3359) 
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VI. TRIPLET ABSORPTION 

In the following we present a detailed description of the calculated triplet absorption 
spectra of octacene, nonacene, and decacene, presented in Figures 8, 9 and 10 of the main 
text. 

1. The first peak always corresponds to the ^ 3 B± excited state of the system, whose wave 
function is dominated by the single excitations \H — ¥ L+ 1) + c.c, irrespective of the 
oligoacene in question, or the Coulomb parameters employed. 

2. The second peak corresponds to a y-polarized transition to the 1 3 A~ excited state, 
for all the oligoacenes, irrespective of the Coulomb parameters employed. The most 
important configuration contributing to the many-particle wave function of the 1 3 A~ 
state, is the double excitation \H — > L; H — 1 — > L + 1). 



3. The nature of the third peak is dependent upon the Coulomb parameters employed 
in the PPP model. For the standard parameter case, this peak corresponds to a 

li g and 3 3 A g 



Iff 



mixture of x- and y-polarized the transitions to states, 2?B^ and 3 3 A~ for octacene 
and nonacene. For octacene, the single excitations \H — ¥ L + 4) + c.c. for the 2 3 B 1 
state, and \H — 1 — > L + 3) + c.c. for 3 3 A~ state contribute the most to the respective 
wave functions. For nonacene, the single excitations \H — ¥ L + 3) + c.c. for the 2 3 B 
state, and the double excitation \H — ¥ L; H — 1 — ¥ L + 1} for the 3 3 A~ states domi- 
nate the corresponding wave functions. However, for decacene, the peak corresponds 
to an x polarized transition to the state 2 3 B^ , whose wave function is dominated 
by single excitations, with the configurations \H — ¥ L + 3) + c.c. contributing the most. 

For the screened parameter case, for all the oligomers, the third peak is due to an 
rr-polarized transition to the state 3 3 B7 whose wave function derives the maximum 
contribution from the single excitations \H — > L + 3} + c.c. 

4. As far as the fourth peak is concerned, with standard parameters it corresponds to an 
rr-polarized transition to the 3 3 -B^ excited state, for all the oligoacenes. For octacene, 
it happens to be the most intense peak of the spectrum, but for nonacene and decacene, 
it is a shoulder to the most intense peak. Double excitations \H — y L; H — > L+3)+c.c. 
contribute the most to the many-particle wave function of this state for octacene, 
whereas for nonacene and decacene, the single excitations \H — 1 — ¥ L + 2) + c.c. 
dominate the wave function. 

In the screened parameter spectrum, the fourth peak is due to a y-polarized transiton 
to the 3 3 A^ state for octacene and nonacene. For octacene, the single excitations 
\H — > L + 5) + c.c. contribute the most to the many particle wave function of this 
state, while for nonacene, the single excitations \H — > L + 6) + c.c. dominate the 
state. However, for decacene, the fourth peak is the second most intense peak of the 
spectrum, corresponding to an a>polarized transition to the Q 3 B± state, whose wave 
function is dominated by the double excitations \H — ¥ L; H — > L + 5) + c.c. 

5. For the standard parameter case, the fifth peak is due to an re-polarized transition 
to the state 4 3 5 l9 for all the oligoacenes. For the case of octacene, it appears as a 
shoulder of the most intense peak (peak IV), while for nonacene and decacene it is the 
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most intense peak. For octacene the wave function of this state is composed mainly of 
single excitations, with configurations \H — 1 — > L + 2) + c.c. contributing the most. 
For nonacene and decacene, however, the wave function is dominated by the double 
excitations \H — y L; H — y L + 4) + c.c. 

In the screened parameter spectrum also the fifth peak corresponds to an x po- 
larized transition for all the oligomers. For octacene, nonacene, and decacene, the 
states involved are 4 3 B{ ' Q 3 B{ ' and 8 3 B^ , respectively. For octacene and nonacene, 
this peak is the second most intense one of the corresponding spectra, and the most 
important configuration contributing to the many-particle wave function of the states 
are the double excitations \H — y L; H — y L + 4) + c.c. However, for decacene, it is a 
relatively weaker feature, with the double excitations \H — 1 — >■ L; H — y L + 6) + c.c. 
dominating the wave function of the state. 

6. The sixth peak in the standard parameter spectrum, is formed by a y-polarized tran- 
sition to the state 9 3 A~ for the case of octacene, whose wave function is dominated 
by the single excitations \H — y L + 5) + c.c. and \H — 3 — y L + 1) + c.c. However, 
for nonacene and decacene, it corresponds to an re-polarized transition to the state 
6 3 B7, whose wave function receives maximum contribution from the triple excitations 
\H->L + l',H-*L + l',H-l->L)+c.c. 

With the screened parameters, this peak is formed by an rr-polarized transition to 
the state 7 3 B± for octacene, with double excitations \H — y L + 1; H — 5 —y L) + c.c. 
dominating its wave function. For nonacene and decacene, the peak is due to mixed x 
and y polarized transitions. For nonacene, states 8 3 B^ and 10 3 A~ form this peak, with 
their wave functions dominated by double excitations \H — > L+l; H— 6 — > L+l)+c.c, 
and \H — 1 — ¥ L + l; H — ¥ L + 2) +c.c, respectively. For decacene, the states involved 
are ll 3 B7 and 14 3 A~ , with their wave functions dominated by the double excitations 
\H — > L + 1; H — 1 — > L + 5) + c.c, and single excitations \H — 1 — > L + 8) + c.c, 
respectively. 

7. With the the standard parameters, peak VII for octacene corrresponds to an x- 
polarized transition to the state 7 3 B~[ with the double excitations \H — > L; H — 2 — y 
L + 3) + c.c. dominating its wave function. However, for nonacene and decacene this 
peak corresponds to a mixed x- and y-polarized transition. For nonacene, the states in- 
volved are 7 3 B± and ll 3 A g , with the double excitations \H —y L+l; H—5 — y L)+c.c, 
and single excitations \H — y L + 5) + c.c, contributing the most to their respective 
wave functions. For decacene, the states forming the peak are 8 3 B^ and 11 3 A~ , of 
which the former is dominated by the double and triple excitations, while the latter 
consists mainly of the single excitations. 

Screened parameter calculations predict peak VII to have a mixed x and y polar- 
ized character for all the oligomers. For the case of octacene, four states ll 3 B7 , 
12 3 Bi , ll 3 A~ , and 12 3 A~ , form this peak and the double excitations dominate the 
wave function of the first three states, while 12 3 A~ is dominated by the single exci- 
tations. For nonacene, two states ll 3 B± and 12 3 A^ shape the peak with the triple 
excitations (\H — y L; H — y L; H — 1 — y L + 2) + c.c.) dominating the wave function 
of the former, and single excitations (\H — 1 — y L + 8) -\- c.c.) that of the latter. For 
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decacene, three states, 13 3 B 



is- 



14 3 Sf ff , and 15 3 A 9 



contribute to the peak, with the 



double excitations dominating all their wave functions. 



In the standard parameter spectrum, peak VIII has a mixed x and y polarized char- 
acter for octacene, but only a y-polarized character for nonacene and decacene. For 
octacene, the peak is formed by states 8 3 B^ and 12 3 A g , of which the wave function 
of the former is dominated by the triple and double excitations, while that of the 
latter by double and single excitations. For nonacene, the state in question is 13 3 A~ 
while for decacene it is 12 3 A~ and wave functions in both the cases are dominated 
by doubly-excited configurations. 



In the screened parameter calculations, peak VIII exists only for octacene and 
nonacene, and is due to an re-polarized transition. For octacene it is formed by 
the state 13 3 i?£ " whose wave function is dominated by both single and double exci- 
tations. For nonacene, the peak is caused by the state 12 3 .Br, whose wave function 
mainly consists of double excitations. 



Table XII: Excited states contributing to the triplet ab- 
sorption spectrum of octacene computed using the MRSDCI 
method coupled with the standard parameters in the PPP 
model Hamiltonian. The table includes many-particle wave 
functions, excitation energies, dipole matrix elements, and 
relative oscillator strengths (ROS) of various states. The ex- 
citation energies are with respect to the ^Az, ground state, 
while the dipole matrix elements, and the ROS, are with re- 
spect to the 1 3 -E>2u state. Below, i +c.c. 1 indicates that the 
coefficient of charge conjugate of a given configuration has 
the same sign, while i —c.c.^ implies that the two coefficients 
have opposite signs. 



Peak 


State 


E (eV) 


Transition ROS 
Dipole (A) 


Wave Functions 


I 


^B- lg 


3.04 


2.950 0.839 


\H-*L + 1) - c.c.(0.5638) 
\H-l->L + 2) +c.c. (0.1836) 


II 


i 3 A- 


3.89 


0.997 0.123 


\H ->• L; H - 1 -> L + 1 ) (0.8187) 
\H -*• L;H - 2 -> L + 2) (0.1135) 


III 


3 3 A g 
^B- lg 


4.28 
4.35 


0.355 0.017 
0.579 0.046 


\H - 1 -► L + 3 ) - c.c.(0.3698) 
\H ^ L + 5) - c.c.(0.3435) 
\H -»• L + 4) - c.c.(0.5362) 

\H-l-*L + 2) +c.c. (0.1999) 


IV 


^B- lg 


5.06 


2.494 1.000 


\H -»• L;H -»• L + 3 ) - c.c.(0.4566) 
\H-l->L + 2) +c.c.(0.3011) 


V 


^B- lg 


5.26 


1.199 0.240 


|ff-l^L + 2)+ c.c.(0.3688) 
\H ^ L;H ^ L + 3) - c.c.(0.2757) 


VI 


#A- 


5.82 


0.588 0.064 


\H ^ L + 5) - c.c.(0.3967) 
\H-3->L + l) -c.c.(0.3809) 
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VII 


^B- lg 


6.10 


1.055 


0.216 


\H -*• L;H - 2 -> L + 3) - c.c.(0.4326) 
\H-+L + 1;H-1-+L + S) - c.c.(0.1274) 


VIII 


&B- lg 


6.29 


0.864 


0.149 


\H-*L + 1',H-*L + 1',H-1-*L) - c.c.(0.4030) 
\H -)■ L; # - 2 ->• L + 3) - c.c.(0.3504) 




12*A- 


6.29 


0.149 


0.004 


|fl'-l-)-L + l;fl'^L + 2) - c.c.(0.3707) 

|#-3^L + 4) -c.c.(0.2968) 



Table XIII: Excited states contributing to the triplet ab- 
sorption spectrum of octacene computed using the MRSDCI 
method coupled with the screened parameters in the PPP 
model Hamiltonian. The rest of the information is same as 





that in 


XI 


I 






Peak 


State 


E (eV) 


Transition 
Dipole (A) 


ROS 


Wave Functions 


I 


^B- lg 


1.97 


4.639 


1.000 


\H-+L + l) +c.c.(0.5936) 
\H - 1 -> L + 2 ) + c.c.(0.0990) 


II 


l*A~ 


2.92 


0.996 


0.064 


\H -> L;H -l^L + 1) (0.8333) 
\H-l-+L;H-l-+L + 2) - c.c.(0.1090) 


III 


&B- lg 


3.39 


0.465 


0.017 


\H ->• L + 3) +c.c.(0.4280) 
|# - 1 -> L + 2 ) + c.c.(0.3614) 


IV 


3 3 AJ 


3.68 


0.621 


0.034 


\H -)• L + 5) +c.c.(0.4915) 
|# - 1 -> L + 4 ) + c.c.(0.3220) 


V 


^B- lg 


4.01 


2.601 


0.641 


|iJ -»• L;H -4- L + 4) - c.c. (0.5605) 
|fT-»-L + l;fl'-l->L + 4) - c.c.(0.1764) 


VI 


7*B- lg 


4.67 


1.017 


0.114 


\H ^ L + 1;H -5^ L) - c.c.(0.5190) 
\H-+L + l',H-l-+L + 4) - c.c.(0.2712) 


VII 


n 3 ^r g 

12 3 A" 


5.13 

5.15 
5.18 
5.19 


0.151 
0.768 
0.345 
0.148 


0.003 
0.072 
0.015 
0.003 


|IT-l-»-L + l;fl'-»-L + 2) - c.c.(0.3311) 

|IT-l-»-L + 8) - c.c. (0.2700) 

|fl'-»-L + l;fl'-l-»-L + 4) - c.c.(0.4913) 

\H-*L + l]H-5-+L) - c.c.(0.2874) 

\H -> L + 1;H -5-* L) - c.c.(0.4987) 

\H ^ L + 1;H -1^ L + 4) - c.c.(0.2693) 

\H - 1 -»■ L + 8 ) + c.c.(0.3959) 

|#^L;#-l^L + 3) - c.c.(0.2687) 


VIII 


i3 3 ^r 9 


5.39 


0.777 


0.077 


|fT - 1 -+ L + 8 ) + c.c.(0.3959) 
\H -> L;H - 1 ->■ L + 3) - c.c.(0.2687) 



Table XIV: Excited states contributing to the triplet ab- 
sorption spectrum of nonacene computed using the MRSDCI 
method coupled with the standard parameters in the PPP 
model Hamiltonian. The rest of the information is same as 
that in IXlT 
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Peak 


State 


E (eV) 


Transition 
Dipole (A) 


ROS 


Wave Functions 


I 


^B~ lg 


2.56 


3.331 


0.785 


\H-+L + 1) +c.c.(0.5522) 
\H -1^ L + 2) +c.c.(0.1825) 


II 


1*A~ 


3.30 


1.018 


0.095 


\H-+L)H-1->L + 1) (0.8148) 
\H -> L;H -l-> L + 3) c.c.(0.1109) 


III 


VA- 


3.79 
3.89 


0.686 
0.352 


0.049 
0.013 


\H -¥ L + 3) -c.c.(0.5156) 

\H-l->L + 2) +c.c.(0.2173) 

\H -> L;H - 1 ->■ L + l) (0.4041) 

|JT-l-»-L + 4) -c.c.(0.3279) 


IV 


^B- lg 


4.39 


0.790 


0.076 


|JT - 1 -> L + 2 ) + c.c.(0.4226) 
\H -> L; H -)• L; iJ - 1 ->• L + 2) + c.c.(0.2222) 


V 


^B- lg 


4.66 


2.789 


1.000 


|if^L;Lf^L + 4) + c.c.(0.5047) 
|F-l^L;#^L + 5) - c.c.(0.2034) 


VI 


&B- lg 


5.14 


0.418 


0.025 


|iT-+L + l;#->-L + l;F-l->-L) +c.c.(0.4111) 

\H -> L; # -» L + 4 ) + c.c.(0.2905) 


VII 


^B- lg 


5.39 

5.44 


0.517 
1.438 


0.040 
0.310 


|#^L + 5) -c.c.(0.3654) 

|JT-l-»-LH-4) -c.c.(0.3219) 

|i!->-L + l;iJ-5->-L) - c.c.(0.4435) 

# -> L + 1; H -»• L + 4) - c.c.(0.3271) 


VIII 


l&A- 


5.75 


0.428 


0.029 


|#-l->-L + l;#->-L + 2) - c.c.(0.2992) 
\H ->■ L; L 7 - 1 ->■ L + 3) - c.c.(0.2305) 



Table XV: Excited states contributing to the triplet absorp- 
tion spectrum of nonacene computed using the MRSDCI 
method coupled with the screened parameters in the PPP 
model Hamiltonian. The rest of the information is same as 







that ii 


iXII 






Peak 


State 


E (eV) 


Transition ROS 
Dipole (A) 


Wave Functions 


I 


l A B- lg 


1.86 


5.218 1.000 


\H->L + 1) +c.c.(0.5919) 
\H - 1 -»• L + 2) - c.c.(0.1143) 


II 


1*A~ 


2.61 


1.032 0.055 


\H-*L]H-1-+L + 1) (0.8266) 
\H - 1 -> L;H -> L + 2) + c.c.(0.1122) 


III 


^B- lg 


3.14 


0.711 0.031 


\H ->• L + 3) +c.c.(0.4278) 
\H-l-*L + 2) -c.c.(0.3551) 


IV 


&A- 


3.58 


0.514 0.019 


\H ^ L + 6) +c.c.(0.4657) 
\H - 1 -»■ L + 4) - c.c.(0.3453) 


V 


^B- lg 


4.09 


2.577 0.535 


\H->L;H-+L + 4) + c.c.(0.5088) 

|L 7 -1^L;L 7 ^L + 1;L 7 ^L + 1) + c.c.(0.2302) 


VI 


WA~ 
&B- lg 


4.52 
4.53 


0. 
1. 


324 0.009 
434 0.184 


|L 7 -l^L + l;L 7 ^L + 2) - c.c.(0.3581) 

\H ^L;H-l^-L + 3) - c.c.(0.2788) 
\H-*L + 1;H~6-+L + 1) - c.c.(0.5621) 
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\H -2^ L + 3) - c.c.(0.1869) 


VII 


^B- g 


4.87 


0.398 


0.015 


\H -+ L;H -> L;H - 1 -» L + 2) - c.c.(0.4361) 
|#-1^L;F-1^L;#^L + 1)+ c.c.(0.3273) 




12% 


4.92 


0.215 


0.004 


|il - 1 -t L + 8 ) + c.c.(0.4643) 
\H -+ L; # - 1 ->• L + 3) - c.c.(0.2204) 


VIII 


l*^ 


5.14 


0.884 


0.079 


|ff-»-L + l;fl'-6-»-L + l) - c.c.(0.5873) 
\H ^ L;H -1^ L + 2) - c.c.(0.0756) 



Table XVI: Excited states contributing to the triplet ab- 
sorption spectrum of decacene computed using the MRSDCI 
method coupled with the standard parameters in the PPP 
model Hamiltonian. The rest of the information is same as 







that i 


tiXII 








Peak 


State 


E (eV) 


Transition 
Dipole (A) 


ROS 


Wave Functions 


I 


^B- lg 


2.49 


3.615 


1.000 


\H-+L + 1) -c.c.(0.5461) 
\H-l-+L + 2) +c.c.(0.1918) 


II 


1*A~ 


3.16 


1.097 


0.117 


\H -+ L;H - 1 ->■ L + l) (0.8189) 
\H-+ L;H-l-> L + 3) - c.c.(0.1177) 


III 


^B- lg 


3.62 


0.800 


0.071 


\H -> L + 3) -c.c.(0.5228) 
\H-l->L + 2) +c.c.(0.1965) 


IV 


^B- lg 


4.35 


0.793 


0.084 


\H - 1 ->• L + 2) + c.c.(0.4745) 
\H-l->L + 4:) -c.c.(0.1760) 


V 


± A B- lg 


4.62 


2.534 


0.912 


\H -+ L;H -► L + 4 ) - c.c.(0.4009) 
|JET -»■ L + 7) -c.c.(0.2788) 


VI 


&B- lg 


5.22 


1.183 


0.225 


|fT-»-L + l;fl'-)-L + l;fl'-l-»-L) - c.c.(0.3604) 
|il -> L + 1; if - 6 -> X ) - c.c.(0.3314) 


VII 


11%- 

&B- lg 


5.42 

5.44 


0.493 
1.179 


0.040 
0.232 


|#^L + 6) +c.c.(0.3537) 

|JT-l-»-LH-4) -c.c.(0.3086) 

\H ^ L + 1; H - 1 ^ L; H - 1 ^ L) - c.c.(0.3928) 

\H-*L + l;H-6-*L) - c.c.(0.3365) 


VIII 


VPA- 


5.70 


0.456 


0.036 


\H-l->L + l;H->L + 2) - c.c.(0.3473) 
\H-+L + l;H-+L + 3) - c.c.(0.2317) 



Table XVII: Excited states contributing to the triplet ab- 
sorption spectrum of decacene computed using the MRSDCI 
method coupled with the screened parameters in the PPP 
model Hamiltonian. The rest of the information is same as 
that in 



XII 



Peak 



State 



E (eV) Transition ROS 
Dipole (A) 



Wave Functions 
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I 


^B- lg 


1.72 


5.841 


1.000 


\H->L + 1) - c.c.(0.5872) 
\H -l-> L + 2) - c.c.(0.1112) 


II 


l^A- 


2.33 


1.074 


0.046 


\H -)• L;H - 1 -> L + 1 ) (0.8198) 

\H ^ L;H -1^ L + 2) - c.c.(0.0986) 


III 


^B- lg 


2.89 


0.831 


0.034 


\H ->• L + 3) - c.c.(0.4352) 
|# - 1 ->• L + 2 ) - c.c.(0.3443) 


IV 


&B- lg 


3.80 


2.695 


0.470 


\H ^ L;H ^ L + 5) - c.c.(0.5362) 
\H->L + l;H-l->L + 5) - c.c.(0.1922) 


V 


^B- lg 


4.24 


1.266 


0.116 


\H -1^ L;H ^ L + 6) - c.c.(0.4998) 
\H->L + l;H-l-+L + 5) - c.c.(0.3005) 


VI 


^'B- lg 


4.70 


0.791 


0.050 


\H->L + l;H-l->L + 5) - c.c.(0.4721) 
\H->L+l;H-6-+L) - c.c.(0.3044) 




IA*A- 


4.72 


0.135 


0.001 


\H-l->L + 8) - c.c. (0.4705) 
\H -3^ L + 5) - c.c.(0.2190) 


VII 


13 3 ^ 


4.90 


1.001 


0.084 


\H->L + l;H-5->L + l) - c.c.(0.3314) 
\H->L + 13) - c.c.(0.3095) 




14^" 


4.91 


0.638 


0.034 


\H ^ L + 13) -c.c.(0.4517) 
\H -> L + 1;H -5 -> L + l) - c.c.(0.2219) 




15MJ 


4.91 


0.162 


0.002 


\H-l->L + l;H->L + 4:) - c.c.(0.4666) 
\H-l->L + l;H-2->L + 2) (0.4449) 
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